Maple 2025 Questions and Posts

What can cause Physics to be activated when one is...

Asked by: nm 11548

Yesterday at 7:32 PM

0 0

In 2025.2, I found that calling Physics:-Setup causes code of `convert/diff` to change (on the fly) causing error.

I removed from my code the places I was calling Physics:-Setup explicitly, but the error still shows up. Which means there is still places in the code to cause Physics to come into play.

I just installed 2025.2 and have not done anything yet. Just running a test.

Does one need to install Physics package?

First, let me give an example showing that calling Physics:-Setup changes `convert/diff` which ends up causing an internal Maple error

>	interface(version);

>	SupportTools:-Version();

>	Physics:-Version();

>

restart;

Original result. No problem

>	print(`convert/diff`);

$proc (a, { useeval::truefalse := false }) local ee, zD; ee := `PDEtools/usediff`(a); zD := indets(ee, 'And(function, satisfies(proc (f) options operator, arrow; (op(0, f))::`PDEtools/D/structure` end proc))'); if zD = {} then ee else _Env_t_in_use := `tools/get_symbols`(a); if useeval then zD := map(`convert/diff/D_to_eval@diff`, zD, ':-useEval'); zD := map(proc (u) options operator, arrow; lhs(u) = subsop(0 = eval, rhs(u)) end proc, zD); subs(zD, ee) else subs(map(`convert/diff/D_to_eval@diff`, zD), a) end if end if end proc$

>	T:=D(x)(infinity):

>	PDEtools:-difforder(convert(T,diff))

After call Physics:-Setup

>

restart;

>	Physics:-Setup('assumingusesAssume' = false):

>	print(`convert/diff`);

proc () local ans; ans := usepdiff(args); if has([ans], 'D') then ans := Physics:-OriginalValue["convert/diff"](ans, args[2 .. -1]) end if; ans end proc

>	T:=D(x)(infinity):

>	PDEtools:-difforder(convert(T,diff))

Error, (in unknown) invalid input: diff received infinity, which is not valid for its 2nd argument

Download why_code_changes_calling_Physics_maple_2025_2_nov_24_2025.mw

Notice that now calling PDEtools:-difforder(convert(T,diff)) fail because Physics:- was called.

But my code now has no calls to Physics:-, I made sure to comment all those out.

But I still get same error when I run my test program. Not able to find where Physics is "activated" again to cause the code change in `convert/diff` that happens dynamically when Physics is loaded.

It must be happining on different call internally. But have no idea why and where.

This is all on new PC. I have no Maple ini files or anything like this that can cause this.

Will installing Physics package explicity fixes this?

Update

I think I know why this happens. Not all fixes made by Edgardo S. Cheb-Terrab in last Physics update are included in Maple 2025.2 !

See https://mapleprimes.com/questions/240843-Why-PDEtoolsSolve-Fail-On-This-Equation where fix is made in Physics v.1880 which I confirmed in Maple 2025.1.

But now I checked and this problem is back in Maple 2025.2. which means that fix was not merged into Maple 2025.2. How could this happen?

Download fixed_in_2025_1_but_missing_in_2025_2.mw

So I am hoping that figuring how to install Physics update in Maple 2025.2 will fix the current problem with convert and also bring that old fix for PDEtools:-Solve back into Maple 2025.2.

This is something that Maplesoft should have taken care of.

Update

Figured how install Physics update package V 1881 and doing so fixed this problem in Maple 2025.2 !

This tells me that all fixes that were in Physics update package were not merged into Maple 2025.2

This worksheet shows that now code works OK, but need to make sure to install Physics update package in Maple 2025.2

>	interface(version);

>	SupportTools:-Version();

>	Physics:-Version();

>

restart;

Original result. No problem

>	print(`convert/diff`);

$proc (a, { useeval::truefalse := false }) local ee, zD; ee := `PDEtools/usediff`(a); zD := indets(ee, 'And(function, satisfies(proc (f) options operator, arrow; (op(0, f))::`PDEtools/D/structure` end proc))'); if zD = {} then ee else _Env_t_in_use := `tools/get_symbols`(a); if useeval then zD := map(`convert/diff/D_to_eval@diff`, zD, ':-useEval'); zD := map(proc (u) options operator, arrow; lhs(u) = subsop(0 = eval, rhs(u)) end proc, zD); subs(zD, ee) else subs(map(`convert/diff/D_to_eval@diff`, zD), a) end if end if end proc$

>	T:=D(x)(infinity):

>	PDEtools:-difforder(convert(T,diff))

>

After call Physics:-Setup. now works OK since Physics update is installed

>

restart;

>	Physics:-Setup('assumingusesAssume' = false):

>	print(`convert/diff`);

proc () local ans; ans := usepdiff(args); if has([ans], 'D') then ans := Physics:-OriginalValue["convert/diff"](ans, args[2 .. -1]) end if; ans end proc

>	T:=D(x)(infinity):

>	PDEtools:-difforder(convert(T,diff))

>

Download why_code_changes_calling_Physics_maple_2025_2_nov_24_2025_V2.mw

Hidden execution group boundaries...

Asked by: vv 14072

November 24 2025

0 0

The new default in Maple 2025 is to hide the Execution Group Boundaries in Worksheet mode.
In order to make them visible one has to use the menu

View > Show/Hide Contents > Execution Group Boundaries F9

or to press the F9 key.

The question: is there a method to have again visible boundaries by default?

(Hiding the boundaries does not seem to be a good idea because an execution group can have (it is not clear why!) multiple prompts.)

MapleToolbox 2025.2 with Matlab 2025b ...

Asked by: MichalKvasnicka 229

November 19 2025

1 1

Is it possible to install and operate Maple toolbox 2025.2 with MATLAB R2025b on Linux (Ubuntu 24.04) platform?

On Ubuntu 24.04 installation process hangs after location of Maple and Matlab folders at installer window.

Any idea what is wrong?

Of course, I know that actual Maple 2025 system requirements mention support only for MATLAB R2024b, but this is an obsolete version of MATLAB.

And finally, Is there any reliable and effective way how to switch between Maple and MATLAB symbolic engine?

Check for Updates not working...

Asked by: C_R 3537

November 15 2025

1 1

I get this message:

Why is it not working using the menu?

Why do I get back ticks when typesetting with HTML...

Asked by: C_R 3537

November 04 2025

0 4

Why are characters printed with additional back ticks when decimal HTML code is used (see output (3) below). lprint((3)) does not print them.

Looks like an exception. Is there a reason for that? Can it be used for something?

Download back_ticks_with_decimal_code.mw

Symbolic calculus w.r.t. ABSTRACT vector, matrix a...

Asked by: sursumCorda 1279

October 28 2025

1 0

Although several similar problems were asked many years ago (see, e.g., the section “Formal linear algebra” here), there appears to be no new progress so far. It is said that such functionalities exists in the Physics package, but I cannot find any corresponding examples.
In short, can Maple at present calculate these examples in terms of symbolic array constructs completely automatically?

In Maple 2025, why don't all the hotkeys show ...

Asked by: qilianshan 65

October 25 2025

3 1

In Maple 2025.1, why don't all the hotkeys show up when I press the Alt key?

Maple 2025.0 is OK，but update to 2025.1 all the hotkeys don't show up when I press the Alt key.Why?

PDE system in L-shaped cavity...

Asked by: KIRAN SAJJAN 50

October 24 2025

0 0

Hello Maple Community,

I'm working on solving a PDE system for fluid flow in an L-shaped cavity (similar to previous work on H-shaped and square domains ) recently asked by some one here. I've implemented the governing equations in Maple but need help generating contour plots.

I'm attaching:
1. My Maple worksheet l_shape_cavity.mw
2. A PDF with the problem description and equations L_shape_cavity_work.pdf

The main issues I'm facing are:
- Setting up proper boundary conditions for the L-shaped domain
- Generating contour plots for velocity and pressure fields
- Ensuring the solution converges properly

Could you please help me with:
1. Correct implementation of the L-shaped domain geometry
2. Generating proper contour plots
3. Any suggestions for improving the numerical solution

Thank you for your assistance!

Is opening recent documents not working when no wo...

Asked by: C_R 3537

October 19 2025

2 2

With the new ribbon interface (and startup-page disabled) I need to open a new worksheet to make opening files by the recent document menu work.

Can someone reproduce this?

How to change the tab key spacing?...

Asked by: NeillSmith 187

October 13 2025

0 0

I use the tab key to control the formatting of lists.

On my computer the tab key spacing under Maple 2025 is a little more than one inch which is too much. I have been through the Maple Help and cannot find how to change the tab key spacing.

Any help will be greatly appreciated.

My current account email is dead so I can only view answers here.

Solving pde system for H-shaped and square domains...

Asked by: Suryakanth 25

October 10 2025

1 7

Hello everyone,

I am working on reproducing Fig. 12 (Isotherms) and Fig. 13 (Streamlines) from the attached paper, cacity_paper_work.pdf which deals with natural convection of GO–MgO/silicone oil hybrid nanofluid inside square and H-shaped cavities under a periodic magnetic field.

Equations (9–11) in the paper describe the steady-state 2D flow and heat transfer in terms of u(x,y), v(x,y), and θ(x,y). I have already set up the full PDE system, parameters, and boundary conditions (BCs) in Maple.

Now, iusing pdsolve(..., numeric, method=fd),

Solve these equations for both square and H-shaped cavities using the given ICs and BCs, and

Generate Isotherm and Streamline plots (similar to Figs. 12 and 13 in the paper).

I have attached:

The PDF (showing equations and figures), and

My Maple worksheet containing the PDE definitions and parameters.

I would really appreciate guidance on:

Whether dsolve can be used effectively for this PDE system,

How to define the BCs correctly for the H-shaped cavity (especially the partially heated wall, θ=1 for 0.4 < y < 0.6), and

How to extract and plot the temperature (θ) and streamfunction (ψ) fields similar to the paper’s plots.

error_in_cavity_work.mw

Why dsolve do not resolve constant of integration ...

Asked by: nm 11548

October 04 2025

0 3

Given an ode with IC. When solution is explicit, Maple resolves the constant of integration as expected and returns solution with no c__1 in it.

But when asked for implicit solution, also with same IC, it now returns solution with c__1 still there.

Is this by design or a bug? Should not constant of integration be resolved using IC in both cases? If unable to solve for c__1 because solution is implicit, should it then not return solution all?

Does this happen in earlier versions of Maple?

>

restart;

>	interface(version);

>	Physics:-Version();

>	SupportTools:-Version();

>

restart;

>	IC:=D(y)(0)=0; ode:=diff(y(x),x)^2+x*y(x)=0: sol:=dsolve([ode,IC]);

>	sol:=dsolve(ode,'implicit');

y(x) = 0, -x^2/((x^3+9*y(x))*(x^2+3*(-x*y(x))^(1/2)))+3*(-x*y(x))^(1/2)/((x^3+9*y(x))*(x^2+3*(-x*y(x))^(1/2)))-c__1 = 0, x^2/((x^3+9*y(x))*(-x^2+3*(-x*y(x))^(1/2)))+3*(-x*y(x))^(1/2)/((x^3+9*y(x))*(-x^2+3*(-x*y(x))^(1/2)))-c__1 = 0

>	#WHY did not resolve constant of integration here?? sol:=dsolve([ode,IC],'implicit');

x^2/((x^3+9*y(x))*(-x^2+3*(-x*y(x))^(1/2)))+3*(-x*y(x))^(1/2)/((x^3+9*y(x))*(-x^2+3*(-x*y(x))^(1/2)))-c__1 = 0, -x^2/((x^3+9*y(x))*(x^2+3*(-x*y(x))^(1/2)))+3*(-x*y(x))^(1/2)/((x^3+9*y(x))*(x^2+3*(-x*y(x))^(1/2)))-c__1 = 0, y(x) = 0

>

Download why_C_still_in_solution_maple_2025_1_oct_4_2025.mw

Physics:-Setup(assumingusesAssume = true): causes ...

Asked by: nm 11548

October 02 2025

0 1

I found very strange behaviour of Maple 2025.1 on Linux.

Same exact code. Calling timelimit(sol,ode) twice on two different solutions. If I do not add

Physics:-Setup(assumingusesAssume = true)

At the start, then both timelimits finish OK. But once Physics:-Setup(assumingusesAssume = true) is added at the start, the second timelimit hangs.

I waited 2 hrs and nothing happens. Maple just freezes. Can't even stop the server from worksheet by clicking on the red button at lower left corner.

This is using latest Physics. Does anyone know why this happens? It seems due to some memory cache issue?

Make sure to save all your work before trying this just in case you have to kill Maple application.

The strange thing, unable to stop the server by clicking on red button or clicking on RESTART KERNEL icon at top, or even clicking on the debuger icon at lower left corner. Only way was to kill Maple itself from Linux command line.

>	interface(version);

>	Physics:-Version();

>	SupportTools:-Version();

>

restart;

Example 1. Not using Physics:-Setup(assumingusesAssume = true): gives NO hang

>

sol_1:=y(x) = 1/2*(3*tan(RootOf(2*3^(1/2)*ln(2*3^(1/6)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*
sec(_Z)^2-8*3^(1/2))^(1/3))-3^(1/2)*ln(4*3^(1/3)-2*3^(1/6)*(9*tan(_Z)*sec(_Z)^2
+9*3^(1/2)*sec(_Z)^2-8*3^(1/2))^(1/3)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*sec(_Z)^2-\
8*3^(1/2))^(2/3))+12*_C2*3^(1/2)+12*3^(1/2)*x+36*I*_C2+36*I*x+6*_Z))^3*3^(1/2)+
9*tan(RootOf(2*3^(1/2)*ln(2*3^(1/6)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*sec(_Z)^2-8*
3^(1/2))^(1/3))-3^(1/2)*ln(4*3^(1/3)-2*3^(1/6)*(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*
sec(_Z)^2-8*3^(1/2))^(1/3)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*sec(_Z)^2-8*3^(1/2))^
(2/3))+12*_C2*3^(1/2)+12*3^(1/2)*x+36*I*_C2+36*I*x+6*_Z))^2+3*3^(1/2)*tan(
RootOf(2*3^(1/2)*ln((3*3^(1/2)*tan(_Z)*sec(_Z)^2-8+9*sec(_Z)^2)^(1/3)+2)-3^(1/2
)*ln((3*3^(1/2)*tan(_Z)*sec(_Z)^2-8+9*sec(_Z)^2)^(2/3)-2*(3*3^(1/2)*tan(_Z)*sec
(_Z)^2-8+9*sec(_Z)^2)^(1/3)+4)+36*I*_C2+36*I*x+12*_C2*3^(1/2)+12*3^(1/2)*x+6*_Z
))+1)^(1/3)*(I*3^(1/2)-1):
ode:=diff(y(x),x)-y(x)^3 = 8:
timelimit(30,odetest(sol_1,ode));

Error, (in simplify/ln/relations) time expired

>

sol_2:=y(x) = -1/2*(3*tan(RootOf(2*3^(1/2)*ln(2*3^(1/6)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)
*sec(_Z)^2-8*3^(1/2))^(1/3))-3^(1/2)*ln(4*3^(1/3)-2*3^(1/6)*(9*tan(_Z)*sec(_Z)^
2+9*3^(1/2)*sec(_Z)^2-8*3^(1/2))^(1/3)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*sec(_Z)^2
-8*3^(1/2))^(2/3))+12*_C3*3^(1/2)+12*3^(1/2)*x-36*I*_C3-36*I*x+6*_Z))^3*3^(1/2)
+9*tan(RootOf(2*3^(1/2)*ln(2*3^(1/6)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*sec(_Z)^2-8
*3^(1/2))^(1/3))-3^(1/2)*ln(4*3^(1/3)-2*3^(1/6)*(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*
sec(_Z)^2-8*3^(1/2))^(1/3)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*sec(_Z)^2-8*3^(1/2))^
(2/3))+12*_C3*3^(1/2)+12*3^(1/2)*x-36*I*_C3-36*I*x+6*_Z))^2+3*3^(1/2)*tan(
RootOf(2*3^(1/2)*ln((3*3^(1/2)*tan(_Z)*sec(_Z)^2-8+9*sec(_Z)^2)^(1/3)+2)-3^(1/2
)*ln((3*3^(1/2)*tan(_Z)*sec(_Z)^2-8+9*sec(_Z)^2)^(2/3)-2*(3*3^(1/2)*tan(_Z)*sec
(_Z)^2-8+9*sec(_Z)^2)^(1/3)+4)-36*I*_C3-36*I*x+12*_C3*3^(1/2)+12*3^(1/2)*x+6*_Z
))+1)^(1/3)*(1+I*3^(1/2)):
timelimit(30,odetest(sol_2,ode));

Error, (in collect) time expired

>

Example 2. using Physics:-Setup(assumingusesAssume = true): second timelimit always hangs

>

restart;

>	Physics:-Setup(assumingusesAssume = true):

>

sol_1:=y(x) = 1/2*(3*tan(RootOf(2*3^(1/2)*ln(2*3^(1/6)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*
sec(_Z)^2-8*3^(1/2))^(1/3))-3^(1/2)*ln(4*3^(1/3)-2*3^(1/6)*(9*tan(_Z)*sec(_Z)^2
+9*3^(1/2)*sec(_Z)^2-8*3^(1/2))^(1/3)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*sec(_Z)^2-\
8*3^(1/2))^(2/3))+12*_C2*3^(1/2)+12*3^(1/2)*x+36*I*_C2+36*I*x+6*_Z))^3*3^(1/2)+
9*tan(RootOf(2*3^(1/2)*ln(2*3^(1/6)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*sec(_Z)^2-8*
3^(1/2))^(1/3))-3^(1/2)*ln(4*3^(1/3)-2*3^(1/6)*(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*
sec(_Z)^2-8*3^(1/2))^(1/3)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*sec(_Z)^2-8*3^(1/2))^
(2/3))+12*_C2*3^(1/2)+12*3^(1/2)*x+36*I*_C2+36*I*x+6*_Z))^2+3*3^(1/2)*tan(
RootOf(2*3^(1/2)*ln((3*3^(1/2)*tan(_Z)*sec(_Z)^2-8+9*sec(_Z)^2)^(1/3)+2)-3^(1/2
)*ln((3*3^(1/2)*tan(_Z)*sec(_Z)^2-8+9*sec(_Z)^2)^(2/3)-2*(3*3^(1/2)*tan(_Z)*sec
(_Z)^2-8+9*sec(_Z)^2)^(1/3)+4)+36*I*_C2+36*I*x+12*_C2*3^(1/2)+12*3^(1/2)*x+6*_Z
))+1)^(1/3)*(I*3^(1/2)-1):
ode:=diff(y(x),x)-y(x)^3 = 8:
timelimit(30,odetest(sol_1,ode));

Error, (in expand) time expired

>	#this below will now hang

>

sol_2:=y(x) = -1/2*(3*tan(RootOf(2*3^(1/2)*ln(2*3^(1/6)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)
*sec(_Z)^2-8*3^(1/2))^(1/3))-3^(1/2)*ln(4*3^(1/3)-2*3^(1/6)*(9*tan(_Z)*sec(_Z)^
2+9*3^(1/2)*sec(_Z)^2-8*3^(1/2))^(1/3)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*sec(_Z)^2
-8*3^(1/2))^(2/3))+12*_C3*3^(1/2)+12*3^(1/2)*x-36*I*_C3-36*I*x+6*_Z))^3*3^(1/2)
+9*tan(RootOf(2*3^(1/2)*ln(2*3^(1/6)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*sec(_Z)^2-8
*3^(1/2))^(1/3))-3^(1/2)*ln(4*3^(1/3)-2*3^(1/6)*(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*
sec(_Z)^2-8*3^(1/2))^(1/3)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*sec(_Z)^2-8*3^(1/2))^
(2/3))+12*_C3*3^(1/2)+12*3^(1/2)*x-36*I*_C3-36*I*x+6*_Z))^2+3*3^(1/2)*tan(
RootOf(2*3^(1/2)*ln((3*3^(1/2)*tan(_Z)*sec(_Z)^2-8+9*sec(_Z)^2)^(1/3)+2)-3^(1/2
)*ln((3*3^(1/2)*tan(_Z)*sec(_Z)^2-8+9*sec(_Z)^2)^(2/3)-2*(3*3^(1/2)*tan(_Z)*sec
(_Z)^2-8+9*sec(_Z)^2)^(1/3)+4)-36*I*_C3-36*I*x+12*_C3*3^(1/2)+12*3^(1/2)*x+6*_Z
))+1)^(1/3)*(1+I*3^(1/2)):
timelimit(30,odetest(sol_2,ode));